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Abstract. This paper analyses in quantitative detail the effect caused by a 
moving mass on a spherical gravitational wave detector. This applies to situations 
where heavy traffic or similar disturbances happen near the GW antenna. Such 
disturbances result in quadrupole tidal stresses in the antenna mass, and they 
therefore precisely fake a real gravitational signal. The study shows that there 
always are characteristic frequencies, depending on the motion of the external 
masses, at which the fake signals are most intense. It however appears that, 
even at those frequencies, fake signals should be orders of magnitude below the 
sensitivity curve of an optimised detector, in likely realistic situations. 



PACS numbers: PACS numbers: 04.80.Nn, 95.55.Ym, 04.30.Nk 

1. Introduction 

Large mass spherical gravitational wave (GW) detectors |, || constitute sensitive 
systems with a real promise of sighting GW events in the frequency range of 1 kHz, 
even with a rather large bandwidth j|| . In such extremely delicate device as an acoustic 
GW detector is f§ |, §, g, ||, any disturbance, whether internal or external, must be 
considered and, if possible, screened out both in hardware and in software ||, |lC| ]. 
Quite recently, a possibility has arisen to build an underground spherical detector in 
the Spanish Pyrinees , located near a road tunnel. 

Heavy traffic in the tunnel is there almost at all times. It is therefore important 
to assess the effect of the passing vehicles close to the antenna. Purely Newtonian 
action is expected on the sphere which exactly fakes a real GW signal, since it shows 
up as a local tide, with its characteristic quadrupole structure. Hope is that such fake 
signals be weak compared to actual GW signals. Otherwise, even if veto control on 
them could easily be exercised, an almost continuous background of spurious signals 
would contaminate the detector output, eventually concealing a real signal. 

The purpose of this paper is to calculate the effect of traffic on the spherical 
detector by means of a point moving mass model, and to present the results of applying 
the theory to a realistic situation. We shall first consider the effect of a mass moving 
with uniform speed past the detector. As we shall see, this causes a very small impact 
on the detector in any reasonable circumstances, and we thus consider next the effect 
of oscillating parts in the passing vehicle, which can be more important. In section || 

§ To whom correspondence should be addressed. 
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Figure 1. Schematics of a passing mass and a spherical GW detector. Distances 
do not match the scales of any realistic example. 



we pose the general problem, then in sections |3j and [l| we fully consider each case; 
section ^ is finally devoted to comment on the results obtained and their practical 
relevance in a real detector site. 



2. The problem 

The situation is best described by the schematic graphic displayed in figure [I]: a 
spherical GW detector has its centre at the origin of coordinates, and a moving mass m 
travels in a straight line past the detector with a given impact parameter \d\. We shall 
assume that the detector's radius R is much smaller than this parameter, i.e., 

R<L\d\, (1) 

and ask which is the Newtonian gravitational action of the moving mass on the GW 
detector. 

Let then <j){x,t) be the gravitational potential generated by the moving mass at 
point x at time t. We shall of course assume it satisfies Poisson's equation 

V 2 0(a;,i) = -AirGp(x,t) (2) 

where G is Newton's constant, and p(x, t) is the mass density distribution of the 
moving mass. We shall also assume this sufficiently small that it can be considered a 
point mass, i.e., 

p( X) t) = mS^(x-b) (3) 
where (see again figure |l|) 



b = b(t) =d + vt 



(4) 
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if one adopts the (arbitrary) convention that the moving mass crosses the point of 
closest distance to the detector, d, at time t = 0, and that it moves with constant 
velocity v. The gravitational potential is thus 

GWs show up in the detector as local tides, therefore the moving mass will cause 
signal confusion inasmuch as it too produces tides. Under the assumption ([!]) above, 
we can calculate the gravitational potential inside the sphere by the approximate 
expansion 

4>{x, t) = 0(0, t) + di<p(0, t) Xi + \ didj^O, t) x lXj +... (6) 

The gravitational accelerations are correspondingly given by the gradient of the 
potential, i.e., a = — V</>, or 

ai(x, t) = ~d t (f>(0, t) Xi - didj<p(0, t)xj + ... (7) 

The first term in the rhs represents the global pull on the sphere, and cannot 
therefore excite its oscillation modes; these are instead excited by the acceleration 
gradients, or tides, given by the second term. These acceleration gradients can be 
converted to a density of forces (force per unit volume) simply by multiplying them 
by the density of the sphere, g, say, to obtain 

if M (x, t) = -g 3^0(0, t) x 3 = gR™ M (t) x 3 (8) 

We readily find 

MM 3Gm ( bdt) bjjt) 1 \ 

where b(t) = |6(f)|. The interaction of the passing mass with the sphere is described 
by the usual equations of elasticity theory 

Q^2-^ 2 u-(\ + fi)V(V-u) = f(x,t). (10) 

The structure of f MM (x,t) is obviously purely quadrupole, so it can be split up 
much in the same way as a GW excitation[J: 

2 

f MM {x,t)= f (m) {x)C^ m \t) (11) 



m=-2 



where 



fl m \x) = gEf^ Xi , C^\t) = g Ef m) if(t) (12) 



with m = — 2, . . . , 2. The matrices E^ m ^ are given by |IJ] 

I Notation and conventions on sphere's modes and response will follow closely those of reference 
henceforth. 
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The sphere's response to the forces (nTl) is thus determined by the series expansion 



u MM ( X> t) = y ^ 



(14) 



which is formally identical to the response to a GW excitation Our concern is now 
to compare this with the system's response to a real GW. We come to it in the next 
section. 



3. Equivalent signal 

A meaningful comparison between the passing mass fake signal and an actual GW 
signal is best set up in frequency domain, as it is in this form that the detector's 
sensitivity is defined. 

Consider first a point on the sphere's surface, x = Rn, where n is a unit outward 
pointing normal, and calculate its radial displacement u(t), i.e., 

u(t) = n ■ u{Rn,t) . (15) 

Because of the form of the sphere's wavefunctions, the Fourier transform of the 
above is easily seen to be || 

u gw (uj) = ~ jr a n2 A n2 (R) iu 2 ~h{uj) L n2 H , (16) 

71=0 

where L n2 (Lo) is the mode's transfer function 

L n2(") = 2 ^ 9 - 1 — — (17) 

with 7„2 the linewidth of the mode, and h{oj) the Fourier transform of the signal h(t), 
defined by 

h{t) = hij(t) n^j (18) 

% Here, C j (t) is a convolution product 

C 'i? ) W= / C<- m \t') smu n2 (t-t')dt' 
Jo 

between the driving term, C' m '(t), and the antenna mode, sinaj n 2(t). In a real system, the latter 
has an additional damping factor, e -7 " 2 ', which gives a system transfer function with a non-zero 
linewidth, 7„2. 
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Likewise, the radial displacement induced at the same location of the sphere by 
the moving mass is easily inferred from equation (14): 



u MM ( W ) = a n2 A n2 (R) R™ m (l;) n lUj L n2 (u>) 



(19) 



n=0 



We can at this point define an equivalent amplitude associated to the fake moving 
mass signal; the natural definition is j| : 



7 MM 
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We thus find 



~ h MM 
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(20) 



(21) 



It is important to stress here that this expression holds for solid spheres as well 
as for hollow spheres, and indeed for a dual sphere, too. This is simply because the 
mode expansion is identical in the GW signal and in the passing mass signal. 

An order of magnitude estimate of the passing mass equivalent signal is 



accomplished by replacing 
i.e., 



in equation (pll) with 
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(22) 



The Fourier transforms ii^ M (w) are expediently calculated adopting a coordinate 
system where the rc-axis is parallel to the velocity vector v, and the y-axis is in the 
direction of the vector b of figure |l} The result is 

z 2 K 2 {z)~2zK 1 (z) -iz 2 K 1 (z) 

-iz 2 K 1 (z) -z 2 K 2 (z) + zK 1 (z) 

zK x {z) 



pMM/ ,\ _ Cm 



(23) 



where d = \d\, and 

z = uj/loq 



LJ 



with the K 's representing modified Bessel functions |13[: 



r(i/ + i/2) 



r(i/2) 

It is now easily seen that 

Gm 



cos zt 



(l+t 2 )-+V2 



7 MM 
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< 7w5- w ~ 2 /(*) 



where 



f{z) = [z A K 2 (z) + (3 + z 2 ) z 2 X 2 (z) - 3« 3 ^i(«)K- 2 (2)] 1/2 
This function has asymptotic behaviours: 



(24) 

(25) 

(26) 
(27) 
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f(z) ~l+(f) , if z « 1 (28) 

f(z) ~ ^z^e" 2 , if z> 1 , (29) 

and has a rather smooth maximum at z ~ 1, as we see in the plot of figure ^|. 

The salient feature of f(z) as regards our present concern is however its 
exponential roll-off in the high frequency range. Take for example as reasonable figures 
a somewhat fast 10 ton truck travelling at a speed of 40 m/s (~ 140 km/h), with an 
impact parameter d = 80 m; this gives a frequency v/d = 0.5 rad/sec, or <~ 0.1 Hz, which 
means we must go up to z ~ 10 4 for frequencies in the kHz range, where a spherical 
detector will be sensitive. But this produces an utterly meaningless h^^ iv ~ xq -4354 
at (jj/2tt = 1 kHz . . . 

Such number clearly indicates that anything will produce more substantial fake 
signals in the detector, provided it has some vigor in the kHz range. An example 




Figure 3. Schematics of a spherical GW detector and a passing mass which, 
simultaneously, oscillates with amplitude A. Again, distances do not match 
realistic scales. 



could be e.g. internal motions of parts of the vehicle, such as the engine pistons, which 
undergo periodic oscillations of a more appreciable magnitude. We come to this next. 

4. Periodic source 

Consider thus that certain internal parts of the passing vehicle oscillate with (angular) 
frequency a and amplitude A — see figure [j| If we represent with m! the mass of one 
such part then the Newtonian potential at a field position x relative to the detector's 
centre is obviously given by 

Cm 1 

where 

b(t) = (vt, 0, 0) , c(t) = (0, A cos at, 0) (31) 

with the axes convention shown in figure p[ The system response will be given in this 
case by the same formulas of section ||, except that the quadrupole matrix R^j M (u>) 
needs to be replaced with 

3Gm / ( [h{t) + Ci {t)} [6j(t) + c 3 -(t)] 1 , 
"** U |6(t)+c(i)|3 ^ |&(*)+c(f)|* 3^" 



and, in particular, the new equivalent signal will be given by 



2 < 4^- 4 i?*° MM M^ MM H , (33) 



by the same arguments which led to equation (E2| 
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The Fourier transforms of the functions R® MM (t) in equation ( |32| ) are too 
complicated to be performed analytically, but there are certain features which can 
be easily inferred by inspection of the defining formulas. Take for example the 
gravitational potential given by equation (pCl) ; we can recast this in the form 



, s GW 



1 2c(x-b) |c| 2 2 
\x-b\ 2 \x-b\- 



(34) 



and perform a power series expansion of the term in square brackets in the rhs on the 
basis that 

|c|<|sb-6| (35) 

i.e., actually assuming that the amplitude A of oscillations of the mass m' is much 
smaller than the impact parameter \d\ — see figure | The series expansion will thus 
consist in a series of ascending powers of \c\ = A cos at, so it will consequently involve 
all the harmonics of the frequency a/27r, clearly with an amplitude which decreases 



with the harmonic order due to the inequality (35). Clearly, the same harmonic 
structure carries over to the higher derivatives of the potential, and in particular to 
the quadrupole moments i?° MM (i). 

Since we must give up the analytic approach of section we shall now estimate 
the Fourier transforms R® mm (uj) by numerical methods. For this we shall make the 
rather natural choice of using Fast Fourier Transform algorithms ]l4| . 

In this approach, first thing we need is to define an appropriate bandwidth to do 
the analysis. Our standard reference will be a dual sphere GW detector, which is the 
best spherical GW antenna we can think of at present — see || . To accurately match 
that reference we must assess the signal intensities up to a frequency of 3000 Hz, hence 
we need a bandwidth of 6000 Hz, including of course negative frequencies. We shall 
thus approximate Fourier transforms by DFT sums: 

N-l 

i?° MM K) = At R™ M ( m At ) ^ mm ' N - w " = nn Nyq /N , (36) 

where N is the total number of sample points, f^Nyq is the Nyquist, or sampling 
frequency, and At is the time interval between successive samples, i.e., 

At = . (37) 

S £Nyq 



We now use the expansion (|36| ) to determine the equivalent signal (|33|) . The result 
is plotted in figure |i[ where a dual sphere sensitivity curve has been added to assess 
the real effect of two instances of fake signals. 

The graphic very clearly shows how far the effect is from causing a problem for 
gravitational wave physics, even though rather exaggerated data have been assumed 
in both of the instances plotted. These are the figures used in the plot: 

Impact parameter: d = 50 metres 

Horizontal speed: i> = 40m/s (<~ 140 km/h) 

Oscillating mass: m! = 100 kg 

Oscillation amplitude: A = 40 cm 
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and the frequencies of oscillation are 100 Hz (6000 rpm) and 150 Hz (9000 rpm) in 
each of the examples, respectively. The results shown in the graph scale with the 
tabulated values as 

- m'A 

^cquiv , 

vd i 

while it is also seen that narrow peaks happen at the harmonic frequencies of 100 Hz 
or 150 Hz — actually of a/2n : as discussed above. 



5. Concluding remarks 

The calculations in this paper show that there is no practical possibility that noise 
generated by nearby traffic, which exactly fakes quadrupole GW noise, disturbs GW 
astronomy in the sensitivity band of the spherical detector, even if this is chosen as 
the best possible dual sphere system. 

While this is particularly clear in the presented plots, we must still stress 
that those results actually constitute upper bounds on the fake signals, due to 
inequalities (|2^) and (|33|). In addition, emphasis should also be put on the fact that 
most of the background signal we see in figure ^ is due to aliasing of the out-of-band 
frequency components. 




300 500 1000 2000 



frequency [Hz] 



Figure 4. GW equivalent signal generated by an oscillating passing by mass, 
and sensitivity curve of a dual sphere GW detector as considered in reference [w . 
The black solid line corresponds to a mass oscillating at 100 Hz, while the the 
red dash-dotted line corresponds to a faster oscillator of 150 Hz. The plot very 
clearly shows that the signal generated by these sources is orders of magnitude 
away from being detectable. 
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Altogether then, it appears that installation of a very sensitive spherical GW 
detector in a site only a few hundred metres from a heavy traffic road is definitely not 
going to represent a problem for the system performance. 
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